Universal Dynamical Decoupling: Two-Qubit States and Beyond 
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Uhrig's dynamical decoupling pulse sequence has emerged as one universal and highly promising 
approach to decoherence suppression. So far both the theoretical and experimental studies have 
examined single-qubit decoherence only. This work extends Uhrig's universal dynamical decoupling 
from one-qubit to two-qubit systems and even to general multi-level quantum systems. In particular, 
we show that by designing appropriate control Hamiltonians for a two-qubit or a multi-level system, 
Uhrig's pulse sequence can also preserve a generalized quantum coherence measure to the order of 
1 + 0(T N+1 ), with only TV pulses. Our results lead to a very useful scheme for efficiently locking 
two-qubit entangled states. Future important applications of Uhrig's pulse sequence in preserving 
the quantum coherence of multi-level quantum systems can also be anticipated. 

PACS numbers: 03.67.Pp, 03.65.Yz, 07.05.Dz, 33.25.+k 
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I. INTRODUCTION 

Decoherence, i.e., the loss of quantum coherence due 
to system-environment coupling, is a major obstacle for 
a variety of fascinating quantum information tasks. Even 
with the assistance of error corrections, decoherence must 
be suppressed below an acceptable level to realize a useful 
quantum operation. Analogous to refocusing techniques 
in nuclear magnetic resonance (NMR) studies, the dy- 
namical decoupling (DD) approach to decoherence sup- 
pression has attracted tremendous interest. The central 
idea of DD is to use a control pulse sequence to effectively 
decouple a quantum system from its environment. 

During the past years several DD pulse sequences have 
been proposed. The so-called "bang-bang" control has 
proved to be very useful [l|, 0, [|[ with a variety of exten- 
sions. However, it is not optimized for a given period T 
of coherence preservation. The Carr-Purcell-Meiboom- 
Gill (CPMG) sequence from the NMR context can sup- 
press decoherence up to 0(T 3 ) 0]. In an approach called 
"concatenated dynamical decoupling" 0, Q , the decoher- 
ence can be suppressed to the order of 0(T N+1 ) with 
2 N pulses. Remarkably, in considering a single qubit 
subject to decoherence without population relaxation, 
Uhrig's (optimal) dynamical decoupling (UDD) pulse se- 
quence proposed in 2007 can suppress decoherence up to 
0{T N+1 ) with only N pulses 0, 0, S|- In a UDD se- 
quence, the jth control pulse is applied at the time 



Tj =Tsin yjw 



), j i. v. 



(1) 



In most cases UDD outperforms all other known DD con- 
trol sequences, a fact already confirmed in two beautiful 
experiments d, [l(| H3 ■ As a dramatic development in 
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theory, Yang and Liu proved that UDD is universal for 
suppressing single-qubit decoherence [l2j]. That is, for a 
single qubit coupled with an arbitrary bath, UDD works 
regardless of how the qubit is coupled to its bath. 

Given the universality of UDD for suppression of 
single-qubit decoherence, it becomes urgent to examine 
whether UDD is useful for preserving quantum coher- 
ence of two-qubit states. This extension is necessary and 
important because many quantum operations involve at 
least two qubits. Conceptually there is also a big differ- 
ence between single-qubit coherence and two-qubit co- 
herence: preserving the latter often means the storage 
of quantum entanglement. Furthermore, because quan- 
tum entanglement is a nonlocal property and cannot be 
affected by local operations, preserving quantum entan- 
glement between two qubits by a control pulse sequence 
will require the use of nonlocal control Hamiltonians. 

In this work, by exploiting a central result in Yang 
and Liu's universality proof [l2| for UDD in single-qubit 
systems and by adopting a generalized coherence mea- 
sure for two-qubit states, we show that UDD pulse se- 
quence does apply to two-qubit systems, at least for pre- 
serving one pre-determined type of quantum coherence. 
The associated control Hamiltonian is also explicitly con- 
structed. This significant extension from single-qubit to 
two-qubit systems opens up an exciting avenue of dy- 
namical protection of quantum entanglement. Indeed, it 
is now possible to efficiently lock a two-qubit system on 
a desired entangled state, without any knowledge of the 
bath. Encouraged by our results for two-qubit systems, 
we then show that in general, the coherence of an arbi- 
trary M-level quantum system, which is characterized by 
our generalized coherence measure, can also be preserved 
by UDD to the order of 1 + 0(T N+1 ) with only N pulses, 
irrespective of how this system is coupled with its envi- 
ronment. Hence, in principle, an arbitrary (but known) 
quantum state of an m-qubit system with M = 2 m levels 
can be locked by UDD, provided that the required control 
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Hamiltonian can be implemented experimentally. To es- 
tablish an interesting connection with a kicked multi-level 
system recently realized in a cold- atom laboratory [l3[, 
we also explicitly construct the UDD control Hamilto- 
nian for decoherence suppression in three-level quantum 
systems. 

This paper is organized as follows. In Sec. II, we first 
briefly outline an important result proved by Yang and 
Liu hj ; we then present our theory for UDD in two-qubit 
systems, followed by an extension to multi-level quantum 
systems. In Sec. Ill, we present supporting results from 
some simple numerical experiments. Section IV discusses 
the implications of our results and then concludes this 
paper. 



II. UDD THEORY FOR TWO-QUBIT AND 
GENERAL MULTI-LEVEL SYSTEMS 



A. On Yang-Liu's Universality Proof for 
Single-Qubit Systems 

For our later use we first briefly describe one central 
result in Yang and Liu's work [12j | for proving the uni- 
versality of the UDD control sequence applied to single- 
qubit systems. Let C and Z be two time-independent 
Hermitian operators. Define two unitary operator U± 
as follows: 



-i[C±(-l)( JV - 1 >Z](T N -T JV _ 1 



x e 



x e 



-i[CTZ](T 2 -T 1 )-i[C±Z]T 1 



(2) 



Yang and Liu proved that for Tj satisfying Eq. ([I]), we 
must have 



u 



l + 0(T JV+i ) 



(3) 



i.e., the product of (u^^j and differs from unity 

only by the order of 0(T N+1 ) for sufficiently small T. In 
the interaction representation, 



Zj(t) 



= e lCt Ze- iCt 



E 

p=0 



^-[C,[C,...[C,Z]}}, (4) 



pi 



p folds 



hence the above expression for U± can be rewritten in 
the following compact form 



ui N) (T) 



-i/ T ±Fjv(t)Zj(t)dt 



(5) 



where T is the final time, % denotes the time-ordering 
operator, and 



F N (t) = (-iy, for te (T S ,T S+1 ). 



(6) 



As an important observation, we note that though Ref. 
[l2| focused on single-qubit decoherence in a bath, Eq. 
(|3|) was proved therein for arbitrary Hermitian operators 
C and Z. This motivated us to investigate under what 
conditions the unitary evolution operator of a controlled 
two-qubit system plus a bath can assume the same form 
as Eq. ©. 



B. Decoherence Suppression in Two-qubit Systems 

Quantum coherence is often characterized by the mag- 
nitude of the off-diagonal matrix elements of the system 
density operator after tracing over the bath. In single- 
qubit cases, the transverse polarization then measures 
the coherence and the longitudinal polarization measures 
the population difference. Such a perspective is often 
helpful so long as its representation-dependent nature is 
well understood. In two-qubit systems or general multi- 
level systems, the concept of quantum coherence becomes 
more ambiguous because there are many off-diagonal ma- 
trix elements of the system density operator. Clearly 
then, to have a general and convenient coherence measure 
will be important for extending decoherence suppression 
studies beyond single-qubit systems. 

Here we define a generalized polarization operator to 
characterize a certain type of coherence. Specifically, as- 
sociated with an arbitrary pure state |^) of our quantum 
system, we define the following polarization operator, 



(7) 



where I is the identity operator. This polarization oper 
ator has the following properties: 

-.2 



(8) 



where \^ ± ) represents all other possible states of the sys- 
tem that are orthogonal to Hence, if the expectation 
value of "Pim is unity, then the system must be on the 
state | if). In this sense, the expectation value of ~P\y) 
measures how much coherence of the | v I')-type is con- 
tained in a given system. For example, in the single-qubit 
case, V\in) measures the longitudinal coherence if is 
chosen as the spin-up state, but measures the transverse 
coherence along a certain direction if \^>) is chosen as a 
superposition of spin-up and spin-down states. Most im- 
portant of all, as seen in the following, the generalized 
polarization operator "Pt-qt) can directly give the required 
control Hamiltonian in order to preserve the quantum 
coherence thus defined. 

We now consider a two-qubit system interacting with 
an arbitrary bath whose sclf-Hamiltonian is given by 
He — Co- The qubits interact with the environment 



via the interaction Hamiltonian HjE = ^c x j + <J J y c y ^ 
<J 3 z c z> j for j = 1,2, where a 3 x , a y , and a J z are the stan- 
dard Pauli matrices, and c QJ - are bath operators. We 



3 



further assume that the qubit-qubit interaction is given 
by H 12 = Y,k,i={x, y ,z} c M°i°f > where the coefficients c M 
may also depend on arbitrary bath operators. A gen- 
eral total Hamiltonian describing a two-qubit system in 
a bath hence becomes 

H = He + Hie + H2E + H12 

= c + (j\c x ,i + <J l y Cy t i + a\c z , x + a x c x<2 

+ cr x< J z c xz + Oy0 2 x C yx + VyVyCyy + u y a z c yz 

+ a l z a 2 c zx + a\a 2 c zy + (j\(j 2 z c zz . (9) 

For convenience each term in the above total Hamiltonian 
is assumed to be time independent (this assumption will 
be lifted in the end). 

Focusing on the two-qubit subspace, the above total 
Hamiltonian is seen to consist of 16 linearly-independent 
terms that span a natural set of basis operators for all 
possible Hermitian operators acting on the two-qubit sys- 
tem. This set of basis operators can be summarized as 

{^Q}i=l,2,- ,16 = {cfe ® <?l}, (10) 

where cr/c, 07 <G {/, <J Xl u yi a z }, with the orthogonality con- 
dition Tvace{XjXk) = ^8jk- But this choice of basis oper- 
ators is rather arbitrary. We find that this operator basis 
set should be changed to new ones to facilitate operator 
manipulations. In the following we examine the suppres- 
sion of two types of coherence, one is associated with 
non-entangled states and the other is associated with a 
Bell state. 



1. Preserving coherence associated with non- entangled 
states 

Let the four basis states of a two-qubit system be |0) = 
I H), |1) = I U>, 12) - I IT), and |3) = | ||). The 
projector associated with each of the four basis states is 
given by 



|0><0| 


= Po = 




+ ^)(1 


+ °\ 


ll)<ll 


= Pi = 




+ ^)(1 


- a. 


|2)<2| 


= P2 = 


> 




+ a \ 


|3> <3j 


= P 3 = 






- a. 



As a simple example, the quantum coherence to be pro- 
tected here is assumed to be "P|o> = 2|0){0| — /. 

We now switch to the following new set of 16 basis 



operators, 

Y x = V {0) =2P -I 

Y 2 = Po + P^^I + al), 

Y, = P0-P1+ 2P 2 = i(/ - al + 2a 2 z ), 
Y 4 = P -P 1 -P 2 +3P 3 

Y 5 = |l)(3| + |3)(l| = i(^-<7^ )) 

Y 6 = ^\l){i\ +l \2,){l\ = \{al - ala% 

Y 7 = \2)(3\ + \Z)(2\ = ±(crl-ala 2 x ), 

Y 8 = -i\2){3\+i\3){2\ = ±(o*-*lo*), 

Y 9 = |1)(2| + |2)(1| = 1(^ + ^), 

Y w = -i|l)(2|+*|2)<l| = i(^-cri«72), 

Y n = |0)(l| + |l)(0| = i(^ + <7>^), 

Y 12 = -i\0){l\+i\l){0\ = ±(o*+*lo*), 

Y 13 = |0)(2| + |2)(0| = i(^ +( r> z 2 ), 

Y u = - t \0)(2\+i\2)(Q\ = ±(o- y +a y a 2 z ), 

Y 15 = |0)<3| + |3><0| = ~(^-^), 

Yi 6 = -i\0)(3\+i\3)(0\ = ±(a 1 x a 2 y +*y x ). (12) 

Using this new set of basis operators for a two-qubit 
system, the total Hamiltonian becomes a linear combi- 
nation of the Yj = 1 — 16) operators defined above, 
i.e., 

16 

where Wj are the expansion coefficients that can contain 
arbitrary bath operators. The above new set of basis 
operators have the following properties. First, the oper- 
ator Yi in this set are identical with T^o) and hence also 
satisfies the interesting properties described by Eq. ([5]). 
Second, 

[Y^Yx] = 0, for J = 1, 2, - • • ,10; 
{Yj,Yi} + = 0, for 3 = 11,12,- •• ,16, (14) 

where [•] represents the commutator and {•}+ represents 
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an anti-commutator. Third, 

10 16 




E A ^ J I E B i y i 



Vi=ll 



16 

E ^j^' 

10 
E 

3=1 

10 

E^- 

3=1 



(15) 



With these observations, we next split the total uncon- 
trolled Hamiltonian into two terms, i.e., H = Ho + H' , 
where 



H = WyYy + W 2 Y 2 + 



W 10 Y W , 



and 



H' = W n Y u + ■■■ + WieYu. 
Evidently, we have the anti-commuting relation 
{Y u H'}+ = 0, 



(16) 
(17) 
(18) 



an important fact for our proof below. 

Consider now the following control Hamiltonian de- 
scribing a sequence of extended UDD 7r-pulses 



N 



H C = Y^ KS(t - Tj 

3=1 



Ti 



(19) 



After the N control pulses, the unitary evolution operator 
for the whole system of the two qubits plus a bath is given 
by {h = 1 throughout) 



u(T) = e -^ ff o+^'K T - T ^)(-iy 1 ) 



x e 



x e 



-i[H +H'](T N -T N - l 



K-iYi) 



■ [H +H'](T 3 -T 2 



(-iYi) 



x e-^+^^-^H-iYx) 



x e 



-ilHo+H']^ 



(20) 



We can then take advantage of the anti-commuting rela- 
tion of Eq. I|18p to exchange the order between (— iYy) 
and the exponentials in the above equation, leading to 

U(T) = (^ Yl ) N e^ H " + ^ NH '^ T ~ T ^ 
-i[iro+(-i) JV " 1 ff'](TV-r JV _i) 



x e 



x e 



x e 



-i[H +H'](T 3 -T 2 ) 
-•i[Ho-H'](T 2 — Ti) 



-i[Ho+H']T! 
\N ~iH T<- 



X e 

{~iYy) N e~ l "°-3 
{-iYy) N U^ ] {T). 



-i/ T F N (t)H'j(t)dt 



(21) 



Here F/v(<;) is already defined in Eq. ([6]), the second 
equality is obtained by using the interaction represen- 
tation, with H'j(t) = e iHat Hie^ iilat , and the last line 



defines the operator U^\t). Clearly, UY y> is exactly 

in the form of defined in Eqs. ((2|) and (O, with 

Ho replacing C and H' replacing Z. This observation 
motivates us to define 



AN) 



U { _ N \T) 



e- lH ° T 3 e 



So -Fs(t)H' t (t)dt 



(22) 



which is completely in parallel with defined in Eq. 

(O. As such, Eq. ([3]) directly leads to 



U { + N) = 1 + 0(T 



(23) 



With Eq. ([23]) obtained we can now evaluate the co- 
herence measure. In particular, for an arbitrary initial 
state given by the density operator pi 1 the expectation 
value of V^o) at time T is given by 

TTace{U(T) Pl U\T)r l0) } 

= Trace{(-zrOM%* («f ) ) t (<n)^|o>} 

= Tr a ce{(-iY 1 ) N uWp i r m ( W W) f (iYy) N } 

= Trace{p^| )} [l + 0(T N+1 )] , (24) 

where we have used V\o) — Yy, Yy = I, and the anti- 
commuting relation between V\o) and H'. Equation (J2H) 
clearly demonstrates that, as a result of the UDD se- 
quence of N pulses, the expectation value of V\q) is pre- 
served to the order of 1 + 0(T N+1 ), for an arbitrary 
initial state. If the initial state is set to be |0), i.e., 
Trace{/»;'P|o)} = 1, then the expectation value of "P\q\ 
remains to be 1 + 0(T N+1 ) at time T, indicating that 
the UDD sequence has locked the system on the state 

|o> = |TT>. 

In our proof of the UDD applicability in preserving the 
coherence "Pim associated with a non-entangled state, the 
first important step is to construct the control operator 
Y\ = and then the control Hamiltonian H c . As is 
clear from Eq. ([5]), each application of the control op- 
erator Yy = V\o) leaves the state |0) intact but induces 
a negative sign for all other two-qubit states. It is in- 
teresting to compare the control operator Yy with what 
can be intuitively expected from early single-qubit UDD 
results. Suppose that the two qubits are unrelated at 
all, then in order to suppress the spin flipping of the first 
qubit (second qubit), we need a control operator <j\ (<xf). 
As such, an intuitive single-qubit-based control Hamilto- 
nian would be 



N 

#c,si„glc = ~ E ~ T 3')( CT 2 + a z 



(25) 



3=1 
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This intuitive control Hamiltonian differs from Eq. (fT9| . 
hinting an importance difference between two-qubit and 
single-qubit cases. Indeed, here the qubit-qubit inter- 
action or the system-environment coupling may directly 
cause a double-flipping error | ft) — > | ||), which can- 
not be suppressed by H c , single- The second key step is 
to split the Hamiltonian H into two parts Hq and H', 
with the former commuting with Y\ and the latter anti- 
commuting with Y±. Once these two steps are achieved, 
the remaining part of our proof becomes straightforward 
by exploiting Eq. 1(25]). These understandings suggest 
that it should be equally possible to preserve the coher- 
ence associated with entangled two-qubit states. 

2. Preserving coherence associated with entangled states 

Consider a different coherence property as defined 
by our generalized polarization operator 'Pim, with \ty) 
taken as a Bell state 

|6) = ^[|U) + |4t)]- (26) 

The other three orthogonal basis states for the two-qubit 
system are now denoted as |1), |2), |3). For example, 
they can be assumed to be |1) = -^=[\ ||) + | J.J.)], 

|2) - ^[| TT) - I 41)], and |3) = ^[| U) - | IT)]- To 
preserve such a new type of coherence, we follow our 
early procedure to first construct a control operator Y\ 
and then a new set of basis operators. In particular, we 
require 

% = P|5> = 2|0><0| - 7 

= \(-i+^>i+^y y -^A)- (27) 

We then construct other 9 basis operators that all com- 
mute with Y\, e.g., 



% = 




% = 


\{I-alal + 2ay y ) : 


% = 


~ a x a x - a y a y ~ 


% = 
% = 

Y 7 = 


l^Pz^x - a x a z)i 

^x - 


Y s = 
% = 




Y 10 = 


1 ^1^2 , „\J2\ 

x y + a y a x)- 



The remaining 6 linearly independent basis operators are 
found to be anti-commuting with Y\. They can be writ- 
ten as 

Yu = \(°i + <%)> 

Y12 = -\{vy z +<?w v )i 

Yl3 = \{°l«l+vWx), 

Y U = ~^l+(T 2 y ), 

Yx 5 = 

Yxe = \{oW y -ay x ). (29) 

The total Hamiltonian can now be rewritten as H = 
H + H' , in which 

Ho = W 1 Y 1 + W 2 Y 2 + ■■■ + W W Y W (30) 

and 

H' = WuYu + ■■■ + W 16 Y 16 . (31) 

It is then evident that if we apply the following control 
Hamiltonian, i.e., 

N y 

H c = ^5(t-T 3 )^ 
i=i 

N 

the time evolution operator of the controlled total system 
becomes entirely parallel to Eqs. ((20]) and pT]) (with an 
arbitrary operator O replaced by O). Hence, using the 
N control pulse described by Eq. f3"2"]) . the quantum 
coherence defined by the expectation value of can 
be preserved up to 1 + 0(T N+1 ), for an arbitrary initial 
state. If the initial state is already the Bell state |0) (i.e., 
coincides with the \^f) that defines our coherence measure 
V\^)), then our UDD control sequence locks the system 
on this Bell state with a fidelity 1 + 0(T N+1 ), no matter 
how the system is coupled to its environment. 

The constant term in the control Hamiltonian H c can 
be dropped because it only induces an overall phase of the 
evolving state. All other terms in H c represent two-body 
and hence nonlocal control. This confirms our initial ex- 
pectation that suppressing the decoherence of entangled 
two-qubit states is more involving than in single-qubit 
cases. 

We have also considered the preservation of another 
Bell state ^-[| U) - 1 IT)]- Following the same procedure 
outlined above, one finds that the required UDD control 
Hamiltonian should be given by 

N 

h c = -X)^(*-W+^+^°2 + °i ff 2X33) 
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which is a pulsed Heisenberg interaction Hamiltonian. 
Such an isotropic control Hamiltonian is consistent with 
the fact the singlet Bell state defining our quantum co- 
herence measure is also isotropic. 



C. UDD in M-level systems 

Our early consideration for two-qubit systems sug- 
gests a general strategy for establishing UDD in an ar- 
bitrary M-level system. Let |0), |1), • • • , \M - 1) be the 
M orthogonal basis states for an M-level system. Their 
associated projectors are defined as Pj = \j)(j\, with 
j = 0, 1, ■ • • , M — 1. Without loss of generality we con- 
sider the quantum coherence to be preserved is of the |0)- 
type, as characterized by V\o) — 2|0)(0| — /. As learned 
from Sec. II-B, the important control operator is then 



Vi = V\ 0) = 2P - I, 



(34) 



with Vi = I. A UDD sequence of this control operator 
can be achieved by the following control Hamiltonian 



H r = 



N 

E 

3=1 



TT5{t 



Tj) — . 

3/ 2 



(35) 



In the M-dimensional Hilbert space, there are totally 
M 2 linearly independent Hermitian operators. We now 
divide the M 2 operators into two groups, one commutes 
with Vi and the other anti-commutes with V\. Specifi- 
cally, the following M — 1 operators 

V 2 = Po + Pi, 

V 3 = P -P 1 + 2P 2 , 



V M 



Po-Pi 



P M _ 2 + (M - l)P M -x (36) 



evidently commutes with V\. In addition, 

other (M — 2)(M — 1) basis operators, denoted 
V M +i,V M +2,- ■ ■ , V M+ (M-2){M-i)> also commute with 
V\. This is the case because we can construct the 
following \{M — 2)(M — 1) basis operators 



\k)(l\ + \l)(k\ 



(37) 



with < k < M and k < I < M. The other 
|(M — 2)(M — 1) basis operators that commute with 
V\ are constructed as 



i\k)(l\ + i\l)(k\, 



(38) 



also with < k < M and k < I < M . All the remaining 
2(M — 1) basis operators are found to anti-commute with 
Vi . Specifically, they can be written as 



Vm+(M-\){M-2)+21-1 
Vm+(M-1)(M-2)+21 

where 1< I < M - 1. 



|0> <Z| + |0(0|; 
-t|0)(l|+i|/)<0|, (39) 



The total Hamiltonian for an uncontrolled M-level sys- 
tem interacting with a bath can now be written as 

Hm = Ho + H , 

M 2 -2M+2 

3=1 
M 2 

h' = E W M> 

j=M 2 -2M+3 



(40) 



where Wj are the expansion coefficients that may contain 
arbitrary bath operators. 

With the UDD control sequence described in Eq. (|3"5|) 
tuned on, the unitary evolution operator can be easily 
investigated using [Ui,7J ] = and {V\,H'} + — 0. In- 
deed, it takes exactly the same form (with Y\ — » Vi) as 
in Eq. (f2"Tj) . We can then conclude that, the quantum 
coherence property Vms associated with an arbitrarily 
pre-selected state \^>) in an M-level system can be pre- 
served with a fidelity 1 + 0(T N+1 ) 1 with only N pulses. 
For an m-qubit system, M = 2"\ In such a multi-qubit 
case, our result here indicates the following: if the initial 
state of an m-qubit system is known, then by (i) setting 
|\&) the same as this initial state, and then (ii) setting 
V\^) as the control operator, the known initial state will 
be efficiently locked by UDD. Certainly, realizing the re- 
quired control Hamiltonian for a multi-qubit system may 
be experimentally challenging. 

Recently, a multi-level system subject to pulsed exter- 
nal fields is experimentally realized in a cold-atom lab- 
oratory [l3j]. To motivate possible experiments of UDD 
using an analogous setup, in the following we consider 
the case of M — 3 in detail. To gain more insights into 
the control operator V\ , here we use angular momentum 
operators in the J — 1 subspace to express all the nine 
basis operators. Specifically, using the eigenstates of the 
j z operator as our representation, we have 



Jx = 



h = IK 





f 1 1 






( ~i 




i 

V * 


f 






(41) 



As an example, we use the state (1,0, 0) T to define our 
coherence measure. The associated control operator V\ 
is then found to be 



vi = j z +r z 



(42) 



Interestingly, this control operator involves a nonlinear 
function of the angular momentum operator j z . This 
requirement can be experimentally fulfilled, because re- 
alizing such kind of operators in a pulsed fashion is one 



7 



main achievement of Ref. [13| , where a "kicked-top" sys- 
tem is realized for the first time. The two different con- 
texts, i.e., UDD by instantaneous pulses and the delta- 
kicked top model for understandin g quan tum-classical 
correspondence and quantum chaos (l3lll4 , can thus 
be connected to each other. 

For the sake of completeness, we also present below 
those operators that commute with V\ , namely, 



V 2 

v 3 
v 5 



1 




+ 2*" 


















+ jzj- 



V2 



{3+3 z -jzj-), 



(43) 



where j± = j x ± ij y ; and those operators that anti- 
commute with Vi, namely, 



V 6 = 

V 7 = 

V a = 

V 9 = 



1 ,. . . \ 

-y=(Jz3++3-3z), 



U-jz - jzj+), 



(44) 



Some linear combinations of these operators will be re- 
quired to construct the control Hamiltonian to preserve 
the coherence associated with other states. 



0.9- 



0.85- 
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— Intuitive pulse, N=8 



'0 0.01 0.02 0.03 



0.04 0.05 OX 
t 



0.07 0.0 



FIG. 1: (color online) Expectation value of the coherence 
measure "Pi*), denoted F(t), as a function of time in dimen- 
sionless units, with being the non-entangled state | ft) 
of a two-qubit system. The bath responsible for the deco- 
herence is modeled by a three-spin system detailed in the 
text. The bottom curve is without any control and the de- 
coherence is significant. The middle curve is calculated from 
a control Hamiltonian intuitively based on two independent 
qubits. The top solid curve represents significant decoherence 
suppression due to our two-qubit UDD control Hamiltonian 
described by Eq. (fT9)l . 



III. SIMPLE NUMERICAL EXPERIMENTS 

To further confirm the UDD control sequences we ex- 
plicitly constructed above, we have performed some sim- 
ple numerical experiments. We first consider a model of 
a two-spin system coupled to a bath of three spins. The 
total Hamiltonian in dimensionless units is hence given 
by 

h = E E h™°? 

m=3 j={x,y,z} 

+E E E E 

n=l k={x,y,z} m>n j={x,y,z} 

+ H c , (45) 

where the first two spins constitute the two-qubit system 
in the absence of any external field, H c represents the 
UDD control Hamiltonian, and the coefficients bj <m and 
Cjk take randomly chosen values in [0, 1] in dimensionless 
units. In addition, to be more realistic, we replace the 
instantaneous S(t — Tj) function in our control Hamil- 

tonians by a Gaussian pulse, i.e., (1 / c^/n)e~^ t ~ Tj " > l c \ 



with c = T/100 unless specified otherwise. Further, we 
set T = 0.1, because this scale is comparable to the de- 
coherence time scale. 

Figure 1 depicts the time dependence of the expecta- 
tion value of the coherence measure "Pi*), denoted F(t), 
with \ty) being the non-entangled state | ft) of the two- 
qubit system. The initial state of the system is also taken 
as the non-entangled state | ft)- As is evident from the 
uncontrolled case (bottom curve) , the decoherence time 
scale without any decoherence suppression is of the or- 
der 0.1 in dimensionless units. Turning on the two-qubit 
UDD control sequence described by Eq. ^TJ]) for N = 8, 
the decoherence (top solid curve) is seen to be greatly 
suppressed. We have also examined the decoherence sup- 
pression using a UDD sequence based on the single-qubit- 
based intuitive control Hamiltonian H c , smg \ c described by 
Eq. <P5)) . As shown in Fig. 1, i? c , single can only produce 
unsatisfactory decoherence suppression. 

Similar results are obtained in Fig. 2, where we aim 
to preserve the coherence measure Vys) associated with 
the Bell state defined in Eq. ([26]) . Apparently, with the 
assistance of our two-qubit UDD control sequence, the 
system is seen to be locked on the Bell state with a fi- 
delity close to unity at all times. Figure 2 also presents 
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FIG. 2: (color online) Same as in Fig. 1, but for V\\n) associ- 
ated with a Bell state defined in Eq. (|26|l . The smooth dashed 
curve represents significant decoherence without control. The 
drastically oscillating dashed curve is calculated from an intu- 
itive single-qubit-based control Hamiltonian, showing strong 
population transfer from the initial state to other two-qubit 
states. The top solid curve represents signficant decoherence 
suppression due to our two-qubit UDD control sequence in 
Eq. (|32}. 



the parallel result if the control Hamiltonian is given by 
He, single shown in Eq. (|25|) . The drastic oscillation of 
F(t) in this case indicates that strong population oscilla- 
tion occurs, thereby demonstrating again the difference 
between single-qubit decoherence suppression and two- 
qubit decoherence suppression. 

Using the same initial state as in Fig. 2, Fig. 3 de- 
picts D = Jq \\p{t) — Pi\\dt, i.e., the time-averaged 
distance between the actual time-evolving density matrix 
from that of a completely locked Bell state, for c = T/100 
and c = T/1000, with different number of UDD pulses. It 
is seen that, at least for the number of UDD pulses con- 
sidered here, c = T/100 = 1/1000 (about one hundredth 
of the decoherence time scale) already suffices to preserve 
a Bell state. That is, there seems to be no need to use 
much shorter pulses such as c = T/1000 = 1/10000, be- 
cause the case of c = T/1000 (dashed line) in Fig. 3 
shows little improvement as compared with the case of 
c = T/100 (solid line). This should be of practical in- 
terest for experimental studies of two-qubit decoherence 
suppression. 

Finally, we show in Fig. 4 the decoherence suppression 
of a three-level quantum system, with the control opera- 
tor given by Eq. (|4"2")) . Here the bath is modeled by other 
four three-level subsystems, and the total Hamiltonian is 



FIG. 3: (color online) The time-averaged distance D between 
the actual density matrix from that of a completely locked 
Bell state, for c = T/100 and c = T/1000, versus the number 
of UDD pulses. The initial state is the same as in Fig. 2. 



chosen as 



H 



m=2 a={x,y,z} 

5 5 

+ Y Y Y Y c <*l3ja,m3p,n 
n=l a={x,y,z} m>n fl={x,y,z} 

+ H c , (46) 



where j a ,m represents the j x , j y , or j z operator associated 
with the mth three-level subsystem, with the first being 
the central system and the other four being the bath. 
The coupling coefficients are again randomly chosen from 
[0, 1] with dimensionless units. The results are analogous 
to those seen in Fig. 1 and Fig. 2, confirming the general 
applicability of our UDD control sequence in multi-level 
quantum systems. Note also that even for the N = 2 
case (middle curve in Fig. 4), decoherence suppression 
already shows up clearly. The results here may motivate 
experimental UDD studies using systems analogous to 
the kicked-top system realized in Ref. 



IV. DISCUSSION AND CONCLUSION 

So far we have assumed that the system-bath cou- 
pling, the bath self-Hamiltonian, and the system Hamil- 
tonian in the absence of the control sequence are all time- 
independent. This assumption can be easily lifted. In- 
deed, as shown in a recent study by Pasini and Uhrig for 
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0.98 
0.96 
0.94 
0.92 

st 0,9 
ST 

0.88 
0.86 
0.84 
0.82 
0.8 



Because the term 3 



-i / T H Q dt 



in Eq. (|4"T)) does not 



■■■N=2 

— N=10 

— Without pulses 
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FIG. 4: (color online) Expectation value of the coherence 
measure denoted F(t), as a function of time in dimen- 

sionless units, with |^) being one basis state of a three-level 
system. The central system is coupled with a bath modeled 
by other four three-level subsystems. The bottom curve rep- 
resents significant decoherence without decoherence control. 
The top two curves represent decoherence suppression based 
on the control operator constructed in Eq. (|42[) . for N — 2 
and N = 10. 



single-qubit systems (l6j, the UDD result holds even af- 
ter introducing a smooth time dependence to these terms. 
The proof in Ref. [l6| is also based on Yang and Liu's 
work [ll| • A similar proof can be done for our extension 
here. Take the two-qubit case with the control operator 
Y\ as an example. If Hq and H' are time-dependent, then 
the unitary evolution operator in Eq. (p?0)) is changed to 



U(T) 



-i ff\H + (-l) N H']dt 



X d 

X 2 

X 2 

X 3 



■ r 1 



z -if£lHa+H']dt 



-i/ Tl [H +H']dt 



-i Jo H dt 



-i^ F N (t)H'j{t)dt 



with 



H'j{t) = J \e l fi Hodt \ H'Z \e-^o H ' 



ndt 



(47) 



(48) 



affect the expectation value of our coherence measure, the 
final expression for the coherence measure is essentially 
the same as before and is hence again given by its initial 
value multiplied by 1 + 0(T N+1 ). 

Our construction of the UDD control sequence is based 
on a pre-determined coherence measure V\q,) that charac- 
terizes a certain type of quantum coherence. This implies 
that our two-qubit UDD relics on which type of deco- 
herence we wish to suppress. Indeed, this is a feature 
shared by Uhrig's work Q and the Yang-Liu universal- 
ity proof [IH for single-qubit systems (i.e., suppressing 
either transverse decoherence or longitudinal population 
relaxation) . Can we also efficiently suppress decoherence 
of different types at the same time, or can we simulta- 
neously preserve the quantum coherence associated with 
entangled states as well as non-entangled states? This is 
a significant issue because the ultimate goal of decoher- 
ence suppression is to suppress the decoherence of a com- 
pletely unknown state and hence to preserve the quan- 
tum coherence of any type at the same time. Fortunately, 
for single-qubit cases: (i) there are already good insights 
into the difference between decoherence suppression for a 
known state and decoherence suppression for an unknown 
state [U (3 (with un-optimized DD schemes); and (ii) a 
very recent study [l9[ showed that suppressing the longi- 
tudinal decoherence and the transverse decoherence of a 
single qubit at the same time in a "near-optimal" fashion 
is possible, by arranging different control Hamiltonians in 
a nested loop structure. Inspired by these studies, we are 
now working on an extended scheme to achieve efficient 
decoherence suppression in two-qubit systems, such that 
two or even more types of coherence properties can be 
preserved. Thanks to our explicit construction of the 
UDD control sequence for non-entangled and entangled 
states, some interesting progress towards this more am- 
bitious goal is being made. For example, we anticipate 
that it is possible to preserve two types of quantum co- 
herence of a two-qubit state at the same time, if we have 
some partial knowledge of the initial state. 

It is well known that decoherence effects on two-qubit 
entanglement can be much different from that on single- 
qubit states. One current important topic is the so-called 
"entanglement sudden death" [2(|, i.e., how two-qubit 
entanglement can completely disappear within a finite 
duration. Since the efficient preservation of two-qubit 
entangled states by UDD is already demonstrated here, 
it becomes certain that the dynamics of entanglement 
death can be strongly affected by applying just very few 
control pulses. In this sense, our results on two-qubit 
systems are not only of great experimental interest to 
quantum entanglement storage, but also of fundamental 
interest to understanding some aspects of entanglement 
dynamics in an environment. 

To conclude, based on a generalized polarization opera- 
tor as a coherence measure, we have shown that UDD also 
applies to two-qubit systems and even to arbitrary multi- 
level quantum systems. The associated control fidelity is 
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still given by 1 + 0(T N+1 ) if N instantaneous control 
pulses are applied. This extension is completely general 
because no assumption on the environment is made. We 
have also explicitly constructed the control Hamiltonian 
for a few examples, including a two-qubit system and a 
three-level system. Our results are expected to advance 
both theoretical and experimental studies of decoherence 
control. 
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